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ABSTRACT 

Direct exoplanet detection is limited by speckle noise in the point spread function (PSF) of the central star. 
This noise can be reduced by subtracting PSF images obtained simultaneously in adjacent narrow spectral 
bands using a multi-channel camera (MCC), but only to a limit imposed by differential optical aberrations in 
the MCC. To alleviate this problem, we suggest the introduction of a holographic diffuser at the focal plane 
of the MCC to convert the PSF image into an incoherent illumination scene that is then re-imaged with the 
MCC. The re-imaging is equivalent to a convolution of the scene with the PSF of each spectral channel of 
the camera. Optical aberrations in the MCC affect only the convolution kernel of each channel and not the 
PSF globally, resulting in better correlated images. We report laboratory measurements with a dual channel 
prototype (1 .575 /im and 1 .625 /.tm) to validate this approach. A speckle noise suppression factor of 12-14 was 
achieved, an improvement by a factor ~5 over that obtained without the holographic diffuser. Simulations of 
realistic exoplanet populations for three representative target samples show that the increase in speckle noise 
attenuation achieved in the laboratory would roughly double the number of planets that could be detected with 
current adaptive optics systems on 8-m telescopes. 

Subject headings: Instrumentation: adaptive optics — planetary systems - 
image processing — techniques: high angular resolution 



stars: imaging — techniques: 



1. INTRODUCTION 

Direct detections of faint exoplanets near bright stars are es- 
sential to extend the census of planetary companions to sepa- 
rations that are beyond the reach of radial velocity surveys and 
to enable measurements of exoplanets' physical properties 
through follow-up multi-color photometry and spectroscopy. 
This difficult endeavour, already being tackled by many, is 
one of the major goals for next generation instruments on 8-m 
class telescopes and future 30- to 100-m telescopes. The main 
difficulty stems from imperfections in the optics that produce 
bright quasi-static speckles in the point spread function (PSF) 
of the central star (Marois et al. 2005; Biller et al. 2004; Mas- 
ciadri et al. 2005). The subtraction of reference PSF images 
is a very efficient way of reducing this speckle noise. A good 
example is angular differential imaging (ADI, Marois et al. 
2006; Lafreniere et al. 2006a), in which images of the target, 
obtained while the field of view rotates, are used for subtract- 
ing the stellar PSF. This technique is among the most success- 
ful at suppressing speckle noise for ground-based imaging; 
however, it is inefficient at small angular separation (< 1") 
because the time required for sufficient natural rotation of the 
field of view is too long. 

Simultaneous spectral differential imaging (SSDI) is a PSF 
subtraction technique that is efficient at all angular separa- 
tions. It consists in the simultaneous acquisition of images in 
adjacent narrow spectral bands within a spectral range where 
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the stellar and planetary spectra differ appreciably (Racine 
et al. 1999; Marois et al. 2000; Smith 1987; Rosenthal et al. 
1996). Judicious image combination and subtraction removes 
the stellar PSF and leaves that of any companion. Note that 
because the PSF images are acquired in different bandpasses, 
they must be rescaled prior to subtraction (diffraction scales 
as A). At larger separations, this scaling allows the technique 
to work even in the absence of differential spectral features 
since a companion would not overlie itself in the rescaled im- 
ages. 

The standard implementation of SSDI is done with a 
multi-channel camera (MCC). The speckle noise suppression 
achievable with such instruments is hampered by differen- 
tial aberrations between the spectral channels (Marois et al. 
2005). For the triple imager TRIDENT (Marois et al. 2005, 
Canada-France-Hawaii Telescope) and the quadruple imager 
SDI (Lenzen et al. 2004, Very Large Telescope), the subtrac- 
tion of two images obtained simultaneously through different 
channels yields a speckle noise attenuation of only ~2-2.5 
(Marois et al. 2005 and Appendix A). Strategies such as mul- 
tiple rotations of the instrument or observation of reference 
stars (Marois et al. 2005; Biller et al. 2004) have been used 
to attenuate further the residuals left by the incomplete si- 
multaneous multi-channel subtraction, and 5er contrast limits 
for magnitude differences of ~9 and > 10 have been achieved 
at an angular separation of 0."5 with TRIDENT and SDI re- 
spectively (Marois et al. 2005; Biller et al. 2006). However, 
the approaches used to circumvent the problem of differential 
aberrations often complicate data acquisition and decrease the 
observing efficiency at the telescope. As several MCCs are 
currently being developed by major observatories for future 
exoplanet searches, e.g. NICI (Toomey & Ftaclas 2003) for 
Gemini South, SPHERE/IRDIS (Dohlen et al. 2006) for the 
VLT, or HiCIAO (Tamura et al. 2006; Hodapp et al. 2006) for 
Subaru, the development of techniques to improve the PSF 
subtraction performance of MCCs is of great interest. 
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In this paper, we propose to minimize the effect of differen- 
tial aberrations by introducing a holographic diffuser (HD) at 
the focal plane of an MCC (see also Lafreniere et al. 2006b). 
This simple approach provides a direct gain in speckle at- 
tenuation without the need to modify the observing strategy. 
The concept and performance estimates are presented in §2. 
Measurements and results using a laboratory prototype are re- 
ported in §3, and §4 presents estimates of the detection limits 
that such a device should reach on 8-m telescopes using cur- 
rent adaptive optics (AO) systems, along with an assessment 
of the corresponding exoplanet detection efficiency. Improve- 
ments in performance for a dedicated camera, implications of 
the laboratory results for other types of SSDI instruments, and 
considerations for implementation of this concept in an exist- 
ing camera are discussed in §5. 

2. MULTI-CHANNEL IMAGER WITH HOLOGRAPHIC 
DIFFUSER 

2.1. Concept 

A light shaping HD is a random surface relief hologram that 
diffuses the light incident upon it into a controlled angular dis- 
tribution. When a coherent wavefront goes through an ideal 
HD, large phase errors are introduced over small scales and 
spatial coherence is lost; thus, every point of the wavefront 
effectively becomes an independent source emitting with an 
angular distribution controlled by the HD. Hence, an HD lo- 
cated at a focal plane converts the PSF image into an inco- 
herent illumination scene, which can then be re-imaged by an 
MCC. This is equivalent to a convolution of the scene with 
the PSF of each channel of the MCC. Optical aberrations in 
the MCC affect only the convolution kernel of each channel 
rather than the PSF globally; the effect of differential aber- 
rations is merely a convolution of the same PSF image with 
different kernels. Hence, contrary to standard MCC imaging, 
the speckle pattern remains the same in all channels and it is 
only the light distribution within each speckle that may vary 
slightly depending on the size of the convolution kernel and 
on the amount of aberrations present in the MCC. If the size 
of the convolution kernel is smaller than the size of a speckle 
in the incident PSF image, then most of the light will remain 
inside the same speckle after convolution and the correlation 
between the different channels will be high. 

An HD is a better choice than other diffuser types for use 
with an MCC because its diffusing cone angle can be made 
small, allowing efficient coupling to relatively large f /# (~8- 
32), typical of existing cameras. Moreover, the diffusing 
properties of an HD are wavelength independent. It should 
be noted that since the diffuser breaks the coherence of the 
wavefront, it precludes the use of a coronagraph downstream 
from it. However, the diffuser poses no problem for a corona- 
graph placed in front of it. 

A real HD does not behave exactly as described above be- 
cause the size of the structures on its surface is not nnegli- 
gible compared to the full-width-at-half-maximum (FWHM) 
of the incident PSF. The surface of an HD can be thought of 
as a random collection of microlenses, each refracting light 
slightly differently such that the overall diffusion of light re- 
sults in the desired angular distribution. A real HD thus trans- 
forms the wavefront incident upon it into a random collection 
of sources, or micro-pupils, of finite size rather than ideal, in- 
finitesimally spaced, point sources. Accordingly, the image 
recorded by an MCC equipped with a focal plane HD would 
correspond to the smooth illumination scene incident upon the 




FIG. 1. — Images of a PSF obtained with a camera equipped with an HD 
at its entrance focal plane. The HD was fixed (left) or moving continuously 
(right) during the exposure; see §2.1 for more detail. The images are shown 
with a logarithmic intensity scale. They were acquired using the experimental 
setup described in §3.1. 

HD broken up into many tiny spots, whose position and inten- 
sity are determined by the position and amount of light inter- 
cepted by the corresponding microlenses (see Figure 1). As 
the micro-pupils produced by the microlenses are usually not 
resolved by the MCC, the shape, or FWHM, of each tiny spot 
is determined by the focal ratio of the MCC. This breakup of 
the smooth illumination scene into a collection of tiny spots is 
problematic, but a simple solution exists. If the HD is moved 
slightly in a plane tangent to its surface, then each tiny spot 
will move on the detector by a corresponding amount, and 
the intensity of each spot will change according to the new 
amount of light intercepted by the corresponding microlens. 
In the limit where the HD moves continuously during an ex- 
posure, the detector is smoothly illuminated and, effectively, 
the image recorded by the MCC corresponds to a convolution 
of the illumination scene incident upon the HD with the PSF 
of the MCC (see Figure 1). 

The imaging process described above can be understood 
mathematically as follows. If I\ and h represent the PSFs of 
the two channels of a dual imager for a point source located 
at its entrance focal plane, and 7; nc represents the PSF inci- 
dent upon the HD, then the resulting image, in each channel, 
is given by the convolution l mc *I\ or Ii nc *h- The residual 
image R obtained by subtracting the images of both channels 
is then given by 

R = (hnc * /l ) - (/inc *h)=Ii BB *(Il-h). ( 1 ) 

The residual image can thus be viewed as being equal to the 
difference of the two MCC PSFs, convolved with the PSF in- 
cident upon the HD. As the spatial scale of the residual noise 
in (I i —I2) is the FWHM of I\ and h, denoted 0mcc> an d the 
width of the convolution kernel is the FWHM of 7j nc , denoted 
f?i nc , the above convolution should reduce the amplitude of 
the residual noise in (I\ - h) by a factor dependent upon the 
ratio ^inc/^MCC- This ratio represents the number of charac- 
teristic lengths over which the residuals are averaged by the 
convolution. A larger ratio 6 mc /6ucc should lead to a larger 
attenuation of the noise. 

It was demonstrated in Marois et al. (2005) that the 
noise attenuation achieved by a dual imager is given by 
ctmcc/Actmcc, where omcc is the amount of aberration in 
each channel and Actmcc is the amount of differential aber- 
ration. Thus when using an HD one can expect, based on 
Eq. (1), that the residual noise also scales as ctmcc/Aomcc* 
but with an additional factor due to the convolution. While 
one should still minimize differential aberrations to maximize 
the attenuation, important leverage is introduced through the 
ratio 6» inc /6>MCC, or equivalently through (//#)i nc /(//#)MCC, 
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TABLE 1 

Speckle noise attenuation from numerical simulations. 



Optical aberrations Att. for S ilK ./S MC c equal to Fit of Eq. (2) 
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80 
80 
80 


50 
50 
50 


1 

2 
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2.1 
4.1 
6.1 


6.9 
13.4 
20.0 


15.2 
29.5 
44.2 


25.3 
49.3 
73.8 


2.1 
2.0 
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1.80 
1.80 
1.81 


80 
80 
80 


25 
25 
25 


1 

2 
3 


4.1 
8.1 
12.2 


16.2 
32.1 
48.0 


37.6 
74.5 
111.6 


65.7 
130.3 
195.2 


4.1 
4.0 
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2.01 
2.01 
2.01 


120 
120 
120 


50 
50 
50 


1 

2 
3 


3.8 
7.4 
11.0 


13.4 
26.0 
38.8 


27.7 
54.1 
81.0 


42.9 
83.7 
125.8 


3.9 
3.8 
3.7 


1.76 
1.76 
1.77 


120 
120 
120 


25 
25 
25 


1 

2 
3 


8.2 
16.2 
24.3 


31.9 
63.0 
94.3 


70.7 
139.9 
209.4 


115.5 
228.2 
342.3 


8.3 
8.2 
8.2 


1.92 
1.92 
1.92 



where (//#)i„ c and (//#)mcc are the focal ratios of the beam 
incident upon the HD and of the MCC. For better perfor- 
mance, these ratios should be made as large as possible. This 
would also minimize the loss in spatial resolution resulting 
from the convolution of /; nc with I\ or h. From convolutions 
of Airy functions, it can be shown that the loss in resolu- 
tion is -12% for (//#W(//#)mcc = 2 and less than 5% for 

(//#)inc/(//#)MCC >3. 

2.2. Speckle noise attenuation estimation 

Numerical simulations were done to estimate the speckle 
noise attenuation that may be achieved with a dual channel 
camera equipped with an HD. For each simulation, three PSFs 
were created: the PSF incident upon the HD (7 inc ) and the 
PSFs of both channels of the MCC (/] and I 2 ). Arrays of 
2048x2048 pixels were used for all calculations. The PSFs 
were obtained as the square modulus of the Fourier Trans- 
form of the complex pupil function. The MCC PSFs FWHM, 
#mcc, was fixed to 5 pixels, and the incident PSF FWHM, # inc , 
was varied to test the effect of the ratio #mc/#MCC- Wavefront 
aberrations having a power-law power spectrum of index -2.7, 
a value appropriate for AO systems and astronomical cameras 
(Marois et al. 2005), were included in the simulations. Differ- 
ent values of the wavefront error of the incident PSF (7i nc , of 
cmcc, and of Aomcc were tried to see their effect. 

Once the three PSFs have been obtained, I\ and h are nor- 
malized to a sum of unity and are subtracted from each other. 
This difference is then convolved with /i nc to obtain the resid- 
ual image. The noise attenuation is finally determined. The 
PSF noise is defined as the standard deviation over an annulus 
of the PSF image after subtraction of an azimuthally averaged 
profile. The noise attenuation r\ is defined as the ratio of the 
PSF noise in one channel over the noise in the residual im- 
age. Many values of #i n c/#MCC, er inc , er M cc, and Aer M cc were 
tried; the results are presented in Table 1 and Figure 2. We 
note that the total amount of static aberration for TRIDENT 
at the CFHT was estimated at 130 nm, of which 50 nm orig- 
inated from the camera, and 120 nm originated from the AO 
system and telescope (Marois et al. 2005); the values of a mc 
and ctmcc tested here are therefore representative of real sys- 
tems. The results indicate that the noise attenuation increases 
with ctmcc / Actmcc and with 0m C /#MCC, in agreement with the 
previous discussion. 

Using a dual channel camera, the maximum value of the 
attenuation, limited by chromaticity of the PSF, is equal to 
A/AA (Marois et al. 2000, 2005), where A is the wavelength 
of the first channel and AA is the difference in wavelength of 
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FIG. 2. — Speckle noise attenuation from numerical simulations as a func- 
tion of inc /8 MC c, the ratio of the FWHM of the PSF incident on the HD 
over that of the intrinsic MCC PSF. The values of <7i nc and <r MCC , in nm, are 
respectively 80 and 50 (crosses), 80 and 25 (stars), 120 and 50 (diamonds), 
and 120 and 25 (triangles). All results are for cr MCC /Acr MCC = 2. The curves 
are the fit of Eq. (2) to each set of points. 

the two channels. The maximum attenuation is thus —30 for 
bandpasses located just outside and just inside the methane 
absorption feature in the spectrum of cold dwarfs at — 1 .6 /zm 
(e.g. Ai — 1 .575/im and \ 2 — 1 .625 /im). From Table 1, the at- 
tenuation should be limited by chromaticity of the PSFs rather 
than by differential aberrations for (//#)i nc /(//#)MCC ^3-4. 

It is found that the attenuation is very well approximated by 
the function 

= a( gMCC ^ ( ^' nc ^ (2) 
VActmcc/ \0mccJ 

where a and (3 are parameters that depend upon <7i nc and (Tmcc 
(see Figure 2); their values are indicated in Table 1. Thus in a 
regime representative of real AO systems and MCC cameras, 
the attenuation that may be achieved with a system equipped 
with an HD roughly scales as the square of 6>i nc / f?MCC- 

3. EXPERIMENTAL RESULTS 

3.1. Dual imager testbed 

A simple testbed operating in the 1.55-1.65 /im wavelength 
range was assembled on an optical bench to validate the con- 
cept presented in the last section. The experimental setup con- 
sists of a PSF generating system followed by a two-channel 
camera. 

The PSF generating imaging system is composed of a 
bi-convex and a plano-convex lens (Thorlabs, LB 1056 and 
LAI 172) and a 6.3 mm diameter pupil stop, which is located 
immediately before the first lens. The system is telecentric 
and has an exit focal ratio of (//#); nc = 64. This focal ratio 
was chosen based on the results of the previous section and the 
// 16 focal ratio of the dual imager described below. The PSF 
is generated by imaging a 25 /im pinhole illuminated from the 
back by a white light source. Due to the small number of op- 
tical components in this system and to the small dimension of 
the beam, the Strehl ratio achieved is very high and no speck- 
les are visible in the resulting PSF image. This is problematic 
for testing a technique aimed specifically at reducing speckle 
noise; it was thus necessary to introduce more aberrations in 
the system. This was achieved by placing an //-band filter 
just before the pupil of the system, resulting in a Strehl ratio 
of —92%, or —75 nm rms of aberrations. 

The exit focal plane of the PSF generating system corre- 
sponds to the entrance focal plane of the two-channel camera. 
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FIG. 3 . — Optical layout of the dual imager testbed. Rays are drawn from two field positions at the entrance focal plane. A separate optical system (not shown) 
is placed in front of the dual imager to produce a PSF on the surface of the HD. See text for more detail. 



An HD can be put in/taken out of the light path at this precise 
location without altering any other part of the setup. For rea- 
sons explained in §2, the HD is mounted inside a ball bearing 
coupled to a small electric motor by a gear belt and spins at 
a rate of ^10 Hz. The rotation axis (center) of the HD is off- 
set by ~6 mm from the optical axis to maximize the effective 
area of the HD seen by the dual imager. The rotation rate was 
chosen to ensure that at least a few complete rotations would 
occur during a single exposure. 

A diagram of the two-channel camera is shown in Figure 3. 
The camera consists primarily of two plano-convex lenses 
(CVI Laser, PLCX-38.1-64.4-C and JML Optical Industries, 
PPX 12925/000); an 8 mm pupil stop is located immediately 
before the second lens. This camera has an entrance focal ra- 
tio of (//#)mcc = 16; this is roughly the smallest focal ratio 
that could be realized easily using simple commercially avail- 
able lenses. A magnification of unity was chosen for sim- 
plicity. Beam splitting is made directly after the second lens 
with a 50/50 beamsplitter cube (Thorlabs, BS012). The beam 
reflected at 90 degrees is folded back toward the detector by 
mean of a right angle prism (Bernhard Halle, UPB0.30); the 
side dimension of this prism was selected such that the ex- 
tra path length in glass compensates precisely the increased 
physical path length of this channel, allowing both channels to 
reach focus in the same plane. The detector used is a Hawaii- 1 
(Rockwell Scientific Company), it is housed inside an IRLab 
cryostat. A square filter mosaic consisting of two abutted rect- 
angular filters of bandwidth 3%, centered on 1.575 /xm and 
1.625 yum respectively, is located in front of the detector and 
allows imaging through different bandpasses in the two chan- 
nels. An additional H-band filter is placed inside the cryostat 
to block radiation longward of 1.8 /im as the narrow band fil- 
ters are not blocked outside of the H band. This dual imager 
has a square field of view 4 mm on a side and a Strehl ratio 
> 95% over the entire field of view in each channel, as esti- 
mated from Zemax. The FWHM of the PSF at the detector 
is 5.5 pixels. There is no anti-reflection coating on any of the 
optical components; hence, many ghosts are present. 

3.2. Data acquisition and reduction 

The speckle noise attenuation capability of the testbed 
MCC was first quantified without the HD. These measure- 
ments were made by acquiring a sequence of images in the 
same bandpass for both channels by rotating the filter mo- 
saic by 90°, which made both channels fall on the same side 
of the mosaic. This eliminates chromatic effects and ensures 
that any decorrelation between the PSFs of the two channels 
is induced by differential aberrations. Then we obtained a se- 



quence of images in the same configuration but with the HD 
in place. This sequence was obtained only a few hours af- 
ter the first one, the only difference being the presence of the 
HD. It thus provides a direct measurement of the improve- 
ment in speckle noise attenuation provided by the HD. Fi- 
nally, a third sequence was acquired with the HD using dif- 
ferent bandpasses for the two channels, which is more repre- 
sentative of astronomical observations. Two diffusers (1° and 
5° FWHM) from Physical Optics Corporation were tried and 
no significant difference in speckle attenuation was found be- 
tween them; however, the 1 diffuser barely fills up the MCC 
pupil stop and yields a transmission of 85% while the 5° one 
largely overfills the MCC pupil stop and yields a transmission 
of 30%. 

All sequences of images were acquired in a similar fashion. 
Five images with the illumination source ON and OFF were 
alternately obtained for a total of 40-50 PSF images and darks. 
Each image consisted of ~5 minutes of coadditions. The in- 
tensity of the illumination source was adjusted to saturate the 
PSF center out to ^3-4 X/D to improve the signal-to-noise 
ratio of the speckles at large separations. 

Basic image reduction consisted in dark subtraction and 
division by the flat field. Bad pixels were replaced by the 
median value of neighboring pixels. For each sequence of 
images, the center of the first image was found by cross- 
correlation with a theoretical PSF having the same FWHM 
in an annulus where the image is not saturated. Then the 
other images of the sequence (including images of the second 
channel) were registered to the first one by cross-correlation 
in the same annulus. The image areas affected by ghost arti- 
facts were masked out for cross-correlation and noise calcu- 
lation purposes. When images were acquired in two band- 
passes simultaneously, the short channel images were spa- 
tially scaled up to compensate for the change in diffraction 
scale with wavelength; the scaling factor was found by cross- 
correlation. In performing the subtraction, the second channel 
PSF intensity was normalized by the ratio of the radial pro- 
files of the two channels. For each sequence, all images and 
all differences were finally median combined. 

3.3. Attenuation results 

A comparison of the images before and after subtraction is 
shown in Figure 4 and the corresponding attenuation curves 
are shown in Figure 5. The single bandpass attenuation with- 
out HD is ^3.5 on average; this is slightly better than, but of 
the order of what is achieved with TRIDENT and SDL The 
addition of the HD boosts the attenuation to ^25, and even 
^40 on Airy rings, an improvement by a factor ^7-10 over 
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FIG. 4. — Comparison of the image of the first channel after subtraction 
of a radial profile (top row) with the difference of the two channels (bottom 
row). The configurations are: single bandpass without HD (left), single band- 
pass with HD (middle), and two-bandpass with HD (right). Display intensity 
range is ±2 X 1CT 5 of PSF peak for all images. The center of the images 
(white) is saturated. Regions affected by the presence of a ghost are masked 
out. 

the configuration without the HD. This is a clear demonstra- 
tion that the effect of differential aberrations is indeed reduced 
by using an HD. Even though the residual noise in this case 
is dominated by pixel-to-pixel noise, some residual speckles 
are present. The attenuation achieved with our testbed may 
be limited by a combination of interpolation errors (when 
shifting the images), flat-field errors, differential pixel cross- 
talk, differential persistence, slight optical distortion of the 
MCC, or differential aberrations. Misalignment of the cube 
and prism may also result in slight field rotation and focus 
differences between the two channels. 

Using the HD and two bandpasses, the speckle noise atten- 
uation achieved is ~12-14. For the two wavelengths of our fil- 
ter mosaic, the maximum attenuation possible for identically 
shaped bandpasses is A/AA ^30. However, the transmission 
profiles of our filters have slopes in opposite directions; this 
decorrelates the two PSFs further. When taking this into con- 
sideration, the maximum attenuation attainable with our mo- 
saic is ~25. Our results show an attenuation lower than this. 
One possible explanation for the difference, in addition to the 
ones mentioned previously, is an inadequate calibration of the 
structures produced by the irregular surface of the HD. Dur- 
ing our measurements, we noticed a slow drift of the detector 
with respect to the pinhole as a function of time. This drift is 
likely correlated with the level of liquid nitrogen in the cryo- 
stat. Hence, the flat field could not be acquired in exactly the 
same configuration as the PSF images and slight calibration 
errors may be present. Such residuals, in the form of arcs 
due to the rotation of the HD, are noticeable in the residual 
image. The calibration of the surface of the HD is important 
only when different bandpasses are used because the images 
have to be rescaled, i.e. a given structure produced by an ir- 
regularity on the HD surface will not spatially coincide with 
itself after rescaling. For this reason, the performance of the 
two-bandpass configuration can be regarded as a lower limit; a 
larger attenuation could be achieved using a more stable setup. 

4. PROJECTED ON-SKY PERFORMANCE 

4. 1 . Detection limit 

Numerical simulations were done to estimate the detec- 
tion limits that currently could be achieved by a dual chan- 
nel camera under three possible scenarios: (1) without an 
HD; (2) with an HD added to an existing MCC - requiring 
the addition of a focal enlarger (§2) that reduces the field of 
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FIG. 5. — Speckle noise attenuation: single bandpass without HD (dotted 
line), single bandpass with HD (dashed line), and two-bandpass with HD 
(solid line). 

view; and (3) with an MCC designed with an HD ab initio - 
avoiding the focal enlarger and the reduction of the field of 
view. For case (2), a 4 x magnification was assumed, yield- 
ing (//#)MCc/(//#)fnc = 4 as in the laboratory experiment. 
The baseline camera was assumed to be similar to NICI, i.e. 
featuring a 1024x1024 pixels detector with a pixel scale of 
0"018 and installed on an 8-m telescope equipped with an 
85-actuator AO system. The simulation was done for a wave- 
length of 1.587 /xm, the central wavelength of the first dif- 
ferential imaging filter of NICI. Coronagraphy was not con- 
sidered as it does not provide significant gain in contrast for 
Strehl ratios typically achieved today, i.e. <50% (Sivaramakr- 
ishnan et al. 2001). 

A pupil of diameter 456 pixels placed at the center of a 
2048x2048 pixels array was used to produce a PSF with a 
FWHM equal to 4.5 pixels; the corresponding pixel scale 
is 0"009 at the wavelength of the simulation. Atmospheric 
phase screens were generated for r ( ) = 20 cm at 0.5 ^m. Addi- 
tionally, a static wavefront phase error having 80 nm rms and a 
power-law power spectrum of index -2.7 was used. AO phase 
filtering was simulated by multiplying the amplitude of the 
phase Fourier transform by a high-pass filter equal to (k / Icao)" 
for k < £ao and to 1 for k > k^o, where k is the spatial fre- 
quency of the phase aberration and k A o is the AO cutoff fre- 
quency (Sivaramakrishnan et al. 2001). The filter index n was 
set to 1.55 to produce an average Strehl ratio of 0.3, the base- 
line value expected for NICI 4 . A long exposure monochro- 
matic PSF was obtained by coadding the result of 5000 in- 
dependent atmospheric wavefront realizations. A polychro- 
matic PSF image was then obtained by coadding 101 spatially 
scaled versions of the monochromatic image, each represent- 
ing a different wavelength within the 1% bandwidth of the 
filter. This PSF image was then either spatially scaled up by 
a factor of 2, or its pixels binned 2 x 2 to obtain proper sam- 
pling for the configuration with the 4 x magnification module 
(0"0045/pixel), or without this module (0."018/pixel), respec- 
tively. The speckle noise profiles were computed from these 
images. 

The residual speckle noise was then obtained by assuming 
attenuation factors of 2.5 and 14 for the configuration without 
and with HD respectively. The pixel-to-pixel noise was com- 
puted for a star of given H-band magnitude assuming a total 
system efficiency of 24% 5 and 19% per channel without and 

4 see the NICI Campaign Science documentation at 
http://www.gemini.edu/sciops/instruments/nici/niciCampaign_orig.html. 

5 This is the estimated total system efficiency of NICI, 
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FIG. 6. — Detection limit in difference of magnitude for an H — 6 star 
and exposure time of 3600 s without HD (solid line), with HD and 4x focal 
enlarger (dashed line), and with HD and no focal enlarger (dot-dashed line). 
The curve for the configuration with the focal enlarger is shorter due to the 
reduced field of view. 

with the HD respectively, a sky background of 14 mag/arcsec 2 
in H, a dark current of 0.15 e~/s per pixel, 30 exposures of 
120 s and a read noise of 10 e"/pixel per exposure. Since 
the statistical distribution of static speckles is non Gaussian, a 
5cr detection limit is not appropriate; rather, to reach a confi- 
dence level similar to that of a 5a Gaussian confidence inter- 
val, the detection limits were obtained as the quadratic sum of 
10(T S peckie and 5(7pi xe i (C. Marois et al., in preparation), where 
""speckle is the speckle noise (dominating at small separations) 
and er pixe i is the pixel-to-pixel noise (dominating at large sep- 
arations). A correction was applied to these limits to account 
for the partial loss in planet flux from the two-channel subtrac- 
tion 6 ; this loss is maximal at zero separation and decreases to 
zero at a separation where the spatial scaling of the images 
displaces sufficiently (by ~1 X/D) the planet. The results for 
an H = 6 star are shown in Figure 6. The higher asymptotic 
limit for the HD with 4x magnification is due to the higher 
relative importance of read noise for this largely oversampled 
configuration. The slightly higher asymptotic limit of the HD 
without magnification compared to the configuration without 
HD is explained by the 20% loss in throughput caused by the 
HD. 

The detection limits presented in Figure 6 should not be re- 
garded as predictions of the best contrast that a dual-channel 
camera could reach in practice since other speckle suppres- 
sion strategies, such as angular differential imaging or sub- 
traction of images obtained after rotation of the instrument, 
could be used to attenuate the noise further. We note however 
that the detection limit without HD presented in this figure is 
very similar to the estimated baseline performance of NICI 7 . 

4.2. Detection efficiency 

A Monte Carlo simulation was done to investigate the im- 
pact of using an HD in terms of actual planet detections. First, 

see http://www.gemini.edu/sciops/instruments/nici/niciCampaign_orig.html. 

6 A maximum loss of 12% was assumed; this is the expected flux loss for a 
T8 dwarf (T e g ~ 800 K) using the narrow band filters of NICI. According to 
the models of Baraffe et al. (2003), giant planets of 5 Mj up or less and older 
than a few 10 Myr should have T cff < 800 K and thus methane absorption 
redward of 1.6 pm is expected to be at least as important for such objects as 
for T8 dwarfs. 

7 see the NICI Campaign Science documentation at 
http://www.gemini.edu/sciops/instruments/nici/niciCampaign_orig.html. 
When comparing both curves, one should keep in mind that the detection 
limits of Figure 6 assume a 10o- S p_._ e criterion whereas the curve in the NICI 
documentation is for 5o- spec kie, this translates into a difference of 0.75 mag. 
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FIG. 7. — Result of the Monte Carlo simulation for the NYS (left), TWA 
(middle), and p Oph (right) samples. The curves show the detection limits 
for a star having an H magnitude equal to the median of the sample for the 
configuration without HD (solid line), with HD and 4x magnification unit 
(dashed line), and with HD and no magnification unit (dot-dashed line). On 
the left and middle panels, the curve for the configuration with the focal en- 
larger is shorter due to the reduced field of view. 

realistic target samples were constructed for three representa- 
tive populations of interest for planet searches. In particular, 
the distance, the age, and the //-band magnitude of each star 
of the samples were determined. The three samples are: 

1. NYS - Nearby young stars (< 150 Myr, < 50 pc). This 
sample includes 72 stars at a distance less than 50 pc 
from the Sun, obtained from Tables 3 and 4 of Wich- 
mann et al. (2003); these stars have ages similar to those 
of the Pleiades stars (20- 150 Myr) based on lithium 
abundance measurements. The stellar distances were 
taken from the same reference. Since the precise age of 
most of these stars is not well known, the stars were ran- 
domly assigned an age from a Gaussian distribution of 
mean 2.0 and standard deviation 0. 15 in log(Age/Myr). 

2. TWA - TW Hydrae Association (~8 Myr, 30-60 pc). 
This sample includes 23 stars from Table 1 in Zucker- 
man & Song (2004). The stellar distances were taken 
from the same reference. The stars were randomly as- 
signed an age from a Gaussian distribution of mean 0.9 
and standard deviation 0.04 in log(Age/Myr). 

3. p Oph - p Ophiuchus star forming cloud (~2 Myr, 
~135 pc). This sample includes 27 stars with R < 15 
from Table 4 of Wilking et al. (2005). The ages were 
taken from the same reference. The stellar distances 
were uniformly sampled in the 120-150 pc interval 
(Knude & Hog 1998). 

In all cases the //-band magnitudes were taken from the 
2MASS point source catalog (PSC, Cutri et al. 2003). For the 
simulation purpose, all samples were artificially increased to a 
total of 50000 targets by re-using the same stars with different 
values for the random parameter. A population of exoplanets 
(one planet per star) was then generated using a mass dis- 
tribution following dN/dM ocM" 1 with 0.5 < M/M Jup < 12, 
a semi-major axis distribution following dN /da oc a" 5 with 
0.1 < a/AU < 50, and a Gaussian eccentricity distribution 
of mean 0.25, standard deviation 0.19 and with < e < 0.8; 
these properties were taken from Marcy et al. (2005) except 
for dN I da which was set according to a minimum mass so- 
lar nebula. All distributions are consistent with the exoplanet 
population from radial velocity surveys. The orbital inclina- 
tion and phase were obtained by uniformly sampling respec- 
tively the orbital orientation and the orbital area swept by the 
planet. Using the stellar distances and the orbital parameters, 
the apparent separation in arcsec was computed for all planets. 
Models of Baraffe et al. (2003) were used to derive the appar- 
ent //-band magnitude of the planets given their mass, age and 
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TABLE 2 
Planet detection efficiency 



Configuration 




Target sample 






NYS 


TWA 


pOph 


No HD 


0.063 


0.175 


0.125 


HD w/4x mag. 


0.113 


0.299 


0.256 


HD w/o mag. 


0.128 


0.295 


0.244 



distance. The magnitude difference between the planets and 
their primary star was finally compared to the detection lim- 
its. The results are shown in Figure 7 and the corresponding 
detection efficiencies are presented in Table 2. The detection 
efficiencies for the configuration without an HD are respec- 
tively 6%, 17.5%, and 12.5% for the NYS, TWA, and p Oph 
samples respectively. The use of an HD roughly doubles these 
numbers. In all cases it is interesting to note that one would 
not suffer from the reduction in FOV needed to use an HD 
with an existing camera, the reason being that most planets 
are located at small angular separations. Given these rela- 
tively small planet detection efficiencies, the relatively small 
number of stars in real target samples, and the unknown but 
likely small fraction of stars with giant planets on large orbits, 
it is reasonable to expect that a dedicated search for planets us- 
ing actual AO systems on 8-m telescopes will uncover at most 
a handful of planets. In this case, i.e. for small number statis- 
tics, the factor of two in detection efficiency resulting from 
the use of an HD could make a significant difference. 

5. DISCUSSION 

The experiment presented in this paper validates the con- 
cept of breaking the wavefront coherence before separating 
it into multiple wavelengths for SSDI. It is clearly demon- 
strated that this strategy reduces the limitations imposed by 
differential aberrations. Still, the two-bandpass attenuation 
obtained with our testbed could have been as high as 25 con- 
sidering the wavelengths of the filter mosaic used. We do not 
believe that the attenuation achieved was limited by differen- 
tial aberrations, nor that it represents a fundamental limit to 
the attenuation achievable with such a system. The fact that 
the single-bandpass attenuation was higher offers evidence for 
this claim. We believe that a dedicated MCC built with more 
care, i.e. with baffles to reduce scattered light, with coat- 
ings on all surfaces to reduce ghost artifacts, with a better 
alignment of the prism with respect to the beam splitter to 
minimize differential rotation and distortion between chan- 
nels, with a more precise control of co-focality between chan- 
nels, and with better mechanical stability, would provide bet- 
ter noise attenuation. The design of NICI correctly addresses 
these issues and thus NICI should not suffer from these short- 
comings. For completeness, we estimated how the planet de- 
tection efficiency would improve with better speckle noise at- 
tenuation. It was previously mentioned that the maximum 
attenuation is ^30 for bandpasses appropriate for exoplanet 
searches. With such an attenuation, the planet detection effi- 
ciency would increase further by ^25%. 

The experimental validation of the concept that we pre- 
sented also offers a good indication that an integral field spec- 
trograph based on a microlens array (e.g. GPI, Macintosh 
et al. 2006a; SPHERE/IFS, Dohlen et al. 2006; or PFI, Mac- 
intosh et al. 2006b) should work well. In such a system, a 
microlens array is placed at the focal plane to spatially sam- 
ple the PSF Each microlens produces a small pupil that is 



dispersed and re-imaged onto a detector by a spectrograph; a 
spectrum of each spatial sample of the PSF is thus obtained. 
Data processing allows the construction of a PSF data cube, 
containing a spectrum for each spatial pixel. Similarly to the 
HD, the spatial sampling of the PSF made by the microlens 
array breaks the coherence of the wavefront and each "micro 
pupil" becomes an independent source for the spectrograph. 
In this case, differential aberrations in the spectrograph af- 
fect only the light distribution within the spectra and not the 
speckles of the PSF imaged. Hence, images of the PSF at 
different wavelengths should be highly correlated. Prelim- 
inary laboratory tests of an integral field spectrograph proto- 
type seem to corroborate this hypothesis (Lavigne et al. 2006). 
Based on a similar argument, the use of a microlens array at 
the entrance focal plane of an MCC, a concept presented in 
Lafreniere et al. (2004), should be a good alternative solution. 

It was previously mentioned that an existing MCC could 
be adapted for use with an HD, we take NICI as an exam- 
ple and briefly state the type of work that would be needed 
to do so. The NICI system (see Toomey & Ftaclas (2003)) 8 
consists of a complete AO system followed by a dual-channel 
camera operating in the 1-5 /im wavelength interval. The fo- 
cal plane of the camera is located outside of the cryostat. At 
this focal plane, a wheel mechanism is used to select between 
various occulting masks for coronagraphy. To implement an 
HD in this system it would thus be necessary to replace the 
occulting mask wheel by a mechanism to hold the HD and 
make it rotate or move continuously as explained in Sect. 2. 
This mechanism would not be technically challenging to re- 
alize as it would be warm. The entrance focal ratio of the 
camera (//#)mcc is 16, thus a diffuser diverging light within 
a cone of half-angle equal to 1.8° would be needed; a custom 
HD matched to this angle should be manufacturable as this 
is in the range of angles available from Physical Optics Cor- 
poration. As explained previously, it would be necessary to 
enlarge the exit focal ratio of the AO system (// 16); based on 
our numerical simulations and experimental results, a factor 
of 3-4 appears adequate. This magnification would be pro- 
vided by a system of lenses placed between the dichroic of 
the AO system and the HD mechanism; these lenses would 
be warm. Accordingly, the field of view would be reduced 
by the same factor. These modifications would preclude the 
use of coronagraphy and would increase the thermal back- 
ground at longer wavelengths. However, the latter is irrel- 
evant for planet searches as SSDI observations are made in 
the //-band, and the loss of coronagraphic capabilities would 
be largely over compensated by the increase in speckle noise 
suppression given the modest Strehl ratios expected. 

6. CONCLUSION 

The concept of placing an HD at the entrance focal plane 
of an MCC to minimize the impact of differential aberrations 
in SSDI was presented. This concept was shown experimen- 
tally to improve significantly the speckle noise attenuation: 
our dual-beam dual-bandpass experimental testbed produced 
a speckle noise attenuation of ~12-14, a factor of ^5 im- 
provement over existing SSDI cameras. Using Monte Carlo 
simulations of realistic exoplanet populations, it was shown 
that such an improvement in speckle noise attenuation could 
yield an increase by a factor ^2 of the number of potentially 
detectable exoplanets. The option presented in this paper 
should be considered seriously for upcoming multi-channel 

8 see also http://www.gemini.edu/sciops/instruments/nici/niciIndex.html 
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cameras and instruments dedicated to exoplanet detection. 
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APPENDIX 

VLT SDI TWO-CHANNEL SUBTRACTION ATTENUATION 

Since the speckle noise attenuation resulting from a two-channel subtraction of images using SDI is unavailable in the literature, 
it is estimated here 9 . Observations of the star AB Dor obtained on 2004 February 1 (Close et al. 2005) were retrieved from the 
European Southern Observatory science archive and used to quantify the performance of SDI. This dataset consists in two sets of 
images at five dither positions obtained at position angles differing by 33 degrees. SDI produces four simultaneous PSF images: 
one at 1.575 /im, one at 1.600 /mi, and two at 1.625 /mi. The reduction steps consisted in the subtraction of a sky frame, division 
by a flat field image and correction of bad pixels using the median of neighboring pixels. The PSF images of all channels were 
then spatially scaled to match the 1.625 /mi channel and registered to a common center with subpixel accuracy. When two PSF 
images were subtracted, their intensity was adjusted to minimize the residual noise. The six possible two-channel subtractions 
were performed and a noise attenuation ^2.-2.5 was obtained for all cases. Figure Al shows the mean attenuation curve over the 
ten images for the subtraction of the two 1.625 /im channels. We stress that this attenuation is not the best achievable with the 
SDI device, subtraction of simultaneous multi-channel difference images obtained at different instrument rotation angles have 
been shown to attenuate further the residuals; see Biller et al. (2004) for more detail. 




o h i i i i H 

10 20 30 40 50 

Angular separation (WD) 

FIG. Al. — Noise attenuation resulting from the subtraction of the two 1.625 (im images with SDI on the VLT (solid line). For completeness, the noise 
attenuation resulting from the subtraction of the 1.58 ^m and 1.625 fim images with TRIDENT at the CFHT is shown as a dashed line; these observations were 
read noise limited beyond 10 X/D (curve adapted from Marois et al. 2005). 
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